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Sections A-B below document some of our recent progress under this project and offer 
several potential routes for the future. 


A. Kinematics of Global Plate Motions and Regional Deformation from VLBI. „ 

Very Long Baseline Interferometry now has the capacity to monitor geodetic positions 
with precisions of a few 1 mm over continental baselines. For tectonic applications, one of 
the major products of the VLBI program is the determination of the rate of change of station 
locations. Vector site velocities are now routinely produced by many groups, including ours. 

We believe that the emphasis in the future will be less on the production vector site velocities 
and more on their interpretation. We are well positioned in this arena, having been in the 
forefront in drawing inferences from VLBI data concerning plate tectonics, plate stability, 
regional deformation and fault slip rates. Below we touch upon one of the novel techniques 
developed by our group in the past two years, VLBI Euler poles. 


A.l. Theory: Euler Poles from VLBI 

Euler poles provide a means to simplify the kinematics of a possibly complicated vector 
field on the surface of a sphere. In the late 1960’s, application of Euler poles to lithospheric 
motions provided a major impetus for plate tectonic theory to sweep the scientific community. 
The use of Euler poles in interpreting space geodetic data has found similar success. 

A typical space geodetic data set produced by NASA (Figure 1) consists of A(« 300) 
length (Rij) and transverse (T,j) velocities determined among A/(~ 30) stations. Reduction 
of the baseline data to a set of vector site velocities is straightforward [Ward, 1988]. Local 
north [t^n,, vjhj] and local cast t^e,] velocities of sites at P, and Pj are linked to the 

observed baseline and transverse T tJ rates by 
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After being weighted by their respective uncertainties, the N equations (1) are inverted to 
obtain 2 M velocity components vj*,vf and their uncertainties cov[uJ*vf]. Rigid plate tectonic 
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The ft, (q = 1, Q) are Euler poles, constant over the m, VLBI sites assigned to the g-th 
plate. 



Fig . L Left . Map illustrating 1 13 VLBI baselines in North America and the Pacific Rim. The dense 
network in the far western United States provides a detailed picture of plate boundary deformation. 

Sites elsewhere in stable North America and the mid-Pacific provide estimates of relative plate 
motion unbiased by plate boundary deformation. 

Fig. 2. Right. Plots of observed V p (P,) (solid arrows) and theoretical V p (P,) (starred vectors) 
site velocities for VM2 at the ten interior sites. VM2 corrects the systematic misfit of NUVEL and 
fits the velocities of the ten stable sites within 2,5 <x. The remaining error reflects the degree of 
consistency of VLBI data and the rigid plate assumption at the ten sites. The triangle, circle, and 
inverted triangle in eastern Canada are the RM2, NUVEL, and VM2 Pacific-North America poles. 

All errors are 3 <r. 

Construction of plate models is an important advancement beyond free-site inversions 
(1). By determining the Euler poles of a few plates upon which several stations ride, bias 
and errors in the velocities of individual sites can be averaged out and failures of the rigid 
plate assumption will be easy to document. In addition, this approach makes it possible to 
quantify the consistency of specified plate models against the VLBI data set as a whole, and 
not just individual baseline rates or site vectors. 

A.2. Global Plate Models. 

One of the original tasks envisioned for space geodesy was the confirmation or confronta- 
tion of geologically- derived plate tectonic portraits with observed instantaneous site motions. 
Progress in this task followed both incremental engineering advances and the slow realization 
of the extent of regional deformation. Through the 1980’s, it became clear that associating a 
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Fig. S. VLBI vector site velocities of European, African and Japanese stations relative to North 
America. Circle labeled 2 in the upper right is the Europe-North American rotation pole from 
NUVEL-1. Inverted triangle labeled 2 is the VLBI pole enclosed in a 3 a error ellipse. Unlike the 
Pacific-North America pole in Figure 2, VLBI and geologic models for European motion disagree. 
Measurements over the next few years should clarify this matter. 



Fig. 4 . Observed site velocities V c (Pi) for western North America and their 3 a error ellipses. 

Lengths of the vectors can be compared with the 20-mm/y standard vector in the box to determine 
site speeds. The observation frame minimizes the velocities of six interior sites on the North America 
plate: Ft. Davis (FD), Flagstaff (F), Platteville (PL), Yuma (Y), Westford, and Fairbanks. The 
stable portion of North America is reached somewhere near the Colorado Plateau. 

site with a lithospheric plate was more complicated than just referencing boundaries drawn 
on a geologic map. For space geodetic data to contrast global tectonic motion models, it is 
imperative to identify stable sites, rigidly attached to the plate in question. 


3 




The Euler pole theory outlined above offers an objective and quantitative means to weed 
out stations which move in a manner inconsistent with the rigid plate assumption. Failure 
of a site^to pass this test may reflect regional or local tectonics or something completely 
unrelated, flexure of an antenna, for instance. Ward [1990] showed that of the 26 United 
States/Pacific sites in the GLB511 solution, only 10 were consistent with rigid two plate 
motion. VLBI-based plate model VM2 (Figure 2), obtained by fitting vector velocities of 
the 10 sites in the stable interiors of the North America and Pacific plates, revealed a misfit 
of VLB I observations and RM2 ( Minster and Jordan , 1978) predictions significant at the 
99% confidence level. NUVEL-1 ( DeMets et al., 1990) reproduced the velocities of these 10 
sites far better than RM2 and offered a fit to the data which is statistically indistinguishable 
from VM2. 

In addition to North America and the Pacific, VLBI is already capable of determining 
plate motions for Europe, Asia and Australia (Figure 3). The observatory mode of VLBI 
is slated to expand through the 1990’s with the North American VLBA and the Soviet 
QUASAR networks. This expansion should permit us to address several tectonic questions of 
global interest including: Asia-North America motion across the arctic, possible Asia-Europe 
expansion, and Asia-Australia/India compression across the Himalayas. In addition, growing 
VLBI ties with South America, Africa and Antarctica should provide further excitement in 
the coming years. 

Our work in global plate motions is not strictly interpretation, but also technique de- 
velopment. Recently we reconsidered the use of VLBI vertical (V) and station height (Z) 
rates of change in equation (2) where historically, only the baseline (B) and transverse (T) 
rates were used to determine local north (N) and (E) velocities. Three new models were 
developed [BT— >NEZ, BTV->NE, BTV-+NEZ] and tested against the BT-+NE standard. 
As reported at the October, 1990 CDP meeting, we found that as yet, the addition of the 
V and Z components fail to provide significant improvements in model consistency. As the 
VLBI vertical component stabilizes over the next few years, we will continue to review its 
status and the procedural underpinnings of our inversion. 

A. 3. Regional Deformation Models. 

As mentioned above, one important application of space geodesy is to identify locations 
where rigid plate tectonics succeeds and where it fails. Judging from the behavior of western 
North America (Figure 4), it seems possible that motions at any continental plate boundary 
may deviate substantially from rigid plate theory. Where stations are sufficiently dense, as in 
western North America, VLBI Euler poles can quantify deformation along specific paths. For 
example, by including velocities of 16 additional stations in the deforming North American 
margin, we have resolved all four terms [Vgjv: Total Pacific-North America motion; Vs^: 
San Andreas slip; V e,w'- Path integrated deformation east and west of the San Andreas] in 
the vector equation characterizing Pacific-North America plate boundary deformation 


Vpn - Vsa = Vg - Vw 

for several points along the San Andreas without reference to a global plate model, geological 
or local geodetic measurements, or restrictions on distributed shear. In most cases, instan- 
taneous VLBI-determined rates agree with geologically-based estimates. VLBI substantially 
eliminates Vw(C), the component of the San Andreas Discrepancy west of C, although a 
few mm/yr of right lateral strike slip on offshore faults can not be rejected. The azimuth 
of San Andreas slip suggests that the transverse ranges north of Los Angeles are absorbing 
8±1 mm/yr of compression perpendicular to the trend of the fault. Convergence across 
the transverse ranges may be slowing, judging from geological studies which infer crustal 
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Fig. 5. Map illustrating the six plates [North America (NAM), Borderland (BLN), Northern 
California (NCA), Southern California (SCA), Vandenberg (VND), Pacific (PAC)] contained in the 
VM models. The North America plate contains the six stations: Ft. Davis, Plattevillc, FlagstafT, 

Yuma, Westford and Fairbanks. The Pacific plate contains the four sites: Vandenberg, Ft. Ord, 
Kwajalein and Kauai. In VM3, Vandenberg and Ft. Ord are detached onto the Vandenberg plate to 
test the need for offshore deformation in central California and west of the Channel Islands. Euler 
poles of the Borderland, Northern California and Southern California plates are used to estimate 
path integrated deformation along NC, N D and N E. 

shortening rates of 17.6 to 26.5 mm/yr over the past two to three million years [Namson and 
Davis , 1988]. Along the Santa Cruz segment of the San Andreas fault which ruptured during 
the 1989 Loma Prieta earthquake, VM2 predicts an azimuth of Pacific- Northern California 
plate motion which is 16 degrees easterly from the local N46°W trend of the fault. This 
discrepancy is compatible with a restraining bend concept for the Santa Cruz segment and 
with the oblique thrust mechanism observed for the quake [Dietz and Ellsworth , 1990]. Even 
a few mm/y of fault normal motion will induce significant rates of vertical uplift and gener- 
ate substantial topography in geologically short order ( Valensise and Ward, 1991). Geodetic 
monitoring of evolving landforms, can, when properly modelled, discriminate modes of tec- 
tonic deformation (i.e. slip on discrete faults versus bulk shortening). In the next several 
years, more precise VLBI height data should play an increasing role. As one example, the 
Palos Verdes mobile VLBI site comes to mind. Palos Verdes is one of the fastest uplifting 
regions in the world, presumably due repeated slip on a buried fault. Continued measure- 
ments on this peninsula over the next several years would be useful in characterizing the 
unseen fault. 

Along plate margins where VLBI sites are not plentiful, path integrated deformation 
vectors are more difficult to obtain' (GPS, of course holds promise here). Nevertheless, it 
still is helpful to classify “normal” deformations from “abnormal” deformations. Ward [1988] 
proposed that the velocities of “normal” sites P,, near the boundary of Plate A and Plate 
B should follow 

V(P.) ={a(P j)Oui + [1 - a{Pi)]n BR }xPi-, 0 < a < 1 (3) 

where D,ar,br are the Euler poles of the plates and a(P) denotes the fractional association 
of P with Plate A. Equation (3) defines a polar velocity field [Minster and Jordan, 1984] 
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Fig. 6 . Lc ft. Cartoon illustrating the behavior of a “normal” plate boundary. For normal sites, the 
VLBI-determined velocity V 0 , should follow equation (3) and fall on the dotted line connecting 
and V b f the predicted rigid plate velocities. Note how boundary deformation can obscure measures 
of relative plate rates. To gauge true relative motion, velocities of sites Pi and P 5 in the stable 
interior must be compared. Including sites P 2 t Psj or P 4 in rigid plate motion estimates would 
induce significant bias. Right. Application of equation (3) to selected plate boundary sites. All of 
these sites are behaving “normally”. 

and simplifies the interpretation of vector deformations to that of characterizing a scalar 
function. Ward [1988] found that velocities (3) fit 88% of the VLBI measurements in the 
western United States when a(P,) was restricted to a two parameter linear transition, and 
95% of the measurements when a(P,) was allowed to be free at each station. 

Figure 6 {left) cartoons a “normal” plate boundary at which a varies monotonically 
across the margin. Points Pi and P5, located in the stable interior of the plates, travel 
at velocities V 0 and V(, predicted by rigid plate theory. Closer to the frontier, rigid plate 
theory fails as the lithosphere deforms. Near the margin, V 0 , the instantaneous site velocity, 
is not likely to be either V Q or V*,, but a mixture described by equation (3). Depending 
on the elastic-plastic composition of the borderlands, strain developing at plate boundaries 
partitions into recoverable and permanent components, if recoverable elastic elements are 
present, then a(Pi) is really a(P, T i), a function of time through the earthquake cycle. 
Averaged over many cycles, the velocities of borderland sites may or may not be equal to 
that observed at any given epoch. Distinctions between geological and instantaneous strain 
accumulation at plate boundaries however, can not affect estimates of Euler poles derived 
from stations in stable plate interiors as these sites have no knowledge of distant earthquake 
cycles. 

Over the next several years we will continue to classify the kinematics of plate margins 
through path integrated deformation vectors where possible (North America, Mediterranean, 
Alaska) or through a(P,) {right, Figure 6). Kinematic classification is the first step toward 
a dynamic analysis. 

B. VLBI and Earthquake Prediction. 

Observations of large scale strain accumulation are indispensable inputs for any master 
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Fig. 7. Individual (Top) and summed ( Bottom ) plate motion vectors V PN , V S a , V E and Vw for 
VM2 projected on the tangent plane at points C, D and E of Figure 5. Error ellipses are 3(7. VLBI- 
determined slip rates at points C,D and E on the San Andreas are 38.6il.3 mm/y (N36 :fc2 W), 

25. Oil. 1 mm/y (N52 # ±3*W), and 10.6i2.1 mm/y (Nl4*iirW) respectively. Segments BLN- 
NAM represent accumulated deformation across the Basin and Range. Segments PAC-BLN, SCA- 
BLN and PAC-NCA represent the San Andreas slip rates at C,D and E. 

model of earthquake prediction. For instance, one would think that the single most important 
number controlling earthquake hazard in the western United States is the rate of relative 
motion between the Pacific and North America Plates. Minster and Jordan’s [1978] RM2 
value of 5.5 cm/yr, for instance, has been built into innumerable models of earthquake 
recurrence. Until recently, rates of plate motion in the historical past could be obtained 
only by indirect methods. Today, space-geodetic techniques not only allow relative plate 
motions to be measured directly, but also reveal where strain is distributed. In Califoinia, 
VLBI has provided the first truly instantaneous estimates of slip rates along various sections 
of the San Andreas (Figure 7) and the first determinations of Euler poles for lithospheric 
sub-plates of the region. In southern California, we find a much slower rate of San Andreas 
slip [25.0 mm/yr (N24°W)] than along the central coast. Most of the decrease [13.3 mm/yr 
(N21°W)] is accommodated in the southern California borderland west of the San Andreas 
to the Channel Islands. Approximately 9 mm/yr of motion taken up by the Basin and 
Range will be transferred to the San Andreas and manifested in steadily increasing slip rates 
south of the Salton Sea. VLBI supports local geodetic evidence ( Prescott et. al . , 1981) for 
a pronounced decrease in San Andreas slip in Northern California. About 15-20 mm/yr of 
Pacific-North America motion is siphoned off to faults east and north of the San Francisco 
Bay. 

VLBI has further suggested several aspects of the tectonic control of earthquake faults. 
When taken together, the six VLBI sites on the Southern California plate (SCA) dictate a 
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Fig. 8. Oblique Mercator projection of California about the VLBI-determined Pacific-Borderland 
pole of rotation (47.18°N, 91.29° W). Because the Borderland has roughly 9 mm /y of N25° W motion 
relative to stable North America, a projection about the PAC-BOR pole is more relevant than the 
more conventional projection about the PAC-NAM pole. The six VLBI sites in southern California 
(stars), when placed on a single rigid plate (SCA), dictate a motion relative to the Borderland 
which is 11° more westerly than PAC-BOR motions, Long and short dashed lines are small circles 
about the PAC-BOR and SCA-BOR poles respectively. 



Fig. 9. Blow up of southern California showing Pleistocene faults. Small circles about the SCA- 
BOR pole (inclined dashed lines) are suggestive of the trends of the Elsinore and San Jacinto faults. 
Loosening of the rigid plate assumption and further localization of relative motions among the 
various faults are priorities for future research. 


8 


■ir 





TIME (YEARS) § STREN 



Fig. 10. Left. Synthetic space-time seismicity ( M > 7) plot for a hypothetical plate boundary 
covering a 900 year time span. The sawtooth curve along the top shows the distribution oflengths 
and strengths for the 18 segments comprising this realization of the boundary. Over about 200 
years, the complex interaction of fault segments and the non-linear failure conditions conspire to 
transform an apparently deterministic model to a chaotic one, in that small changes in initial 
conditions ultimately lead to very different patterns of earthquakes. Right. Plot showing the 
accumulation of tectonic displacement along strike. Each pillow represents the slip resulting from 
one of the discrete events to the left. Averaged over time, all of the segments must maintain a plate 
tectonic rate of 7 cm/yr. Note how some segments (#15, for instance) display a greater fraction of 
characteristic earthquakes ( Schwartz and Coppersmith , 1984) than others (eg. #4). 

sense of motion relative to the North America Borderland (BOR) which is 11° more westerly 
than Pacific-Borderland motion (Figure 8). Most likely in an effort to minimize obliquity 
(angular mismatch of relative fault displacement and fault strike), the Elsinore and San 
Jacinto faults (Figure 9) have aligned themselves along small circles about the SCA-BOR 
pole. On the San Andreas, obliquity migrates from fault-normal tension near the Salton Sea 
to fault-normal compression east of Los Angeles. The details of how the obliquity of SCA- 
BOR and PAC-SCA motions are accommodated remain a major concern for earthquake 
hazard evaluation (eg. blind thrusts). It is a concern in which VLBI will play an important 
role for the foreseeable future. 

The pattern of earthquakes along even a simple plate boundary can be complex, and 
there is growing evidence that the seismic moment release may be only a fraction of the 
total slip budget. Lately, we have made strides in modeling synthetic seismicity and fault 
interaction (Figure 10). These models, based on static dislocation theory ( Bonafcde cl. al. y 
1985; Dragoni , 1990), assist in localizing geodetic experiments and provide guides for the 
magnitude of the signals to be expected. The slip patches in Figure 10 could occur seismically 
or aseismically. VLBI has previously detected coseismic signals associated the the 1989 Loma 
Prieta earthquake ( Clark ct al. y 1990) and the 1988 Alaskan sound event, and should have no 
trouble detecting any M > 7 slip event, even if it occurs aseismically. Combining synthetic 
seismicity models with observed motions of the Alaskan sites, for example, should clarify the 
partitioning of strain transfer among the region’s great strike slip and thrust faults systems 
(Figure 11). Our speculation is that the high velocity of Yakataga is largely permanent, the 
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Fig. 11. Map of VLBI-determined site velocities in Alaska relative to North America. We suspect 
that the high velocity of the Yakataga site reflects permanent strain taken up on interior faults. 

motion being taken up inboard. 

The VLBI network is continuously expanding and improving. The data employed in 
Word [1990] is already two years old. In that interval, there have been several reoccupations 
including: Vernal Utah, Kwajalein and the Alaskan sites. Recovering the motion of Vernal 
is key toward identifying the eastward limit of Basin and Range expansion. Refinement of 
the Kwajalein velocity will, oddly enough (through reduction of the uncertainty in Pacific 
plate stability), help put to rest the question of the slip rate on the San Gregorio- Hosgri 
fault. (Both marine terrace data and VLBI suggest to us that this is probably a fossil 
feature.) With continued NASA funding we will update velocity estimates with the newest 
measurements and attempt to further isolate the deformation field through the introduction 
of multiple or non-rigid tectonic plates. 

C. IMPORTANCE OF MOBILE VLBI IN THE 1990’s. 

The successes we now enjoy in interpreting tectonics from VLBI data have been ten 
years in the making. Much of the progress, certainly in the area of regional deformation, has 
to be credited to the mobile VLBI program. At the March 1991 working group meeting in 
Pasadena, NASA headquarters expressed sentiments that mobile VLBI has outlived its use- 
fulness and indicated that it will not play any role in the CDP follow-on program. I consider 
the elimination of mobile VLBI to be a sad turn of events. The costs saved in terminating 
mobile VLBI have to be weighted against the likely scientific return. In geodetic monitoring 
of tectonics, continuity is critical. The long history and consistency of VLBI observations at 
many sites through the western United States and Alaska will not be duplicated by replace- 
ment technologies for many years, if ever. Mobile VLBI measurements touch upon almost 
every aspect of the work proposed here. If these lines of research are viewed as a valuable 
component of the Dynamics of the Solid Earth NRA, then I would argue strongly that mobile 
VLBI campaigns must continue at some level through the 1990s. 
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